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20.  Abstract 

where  V is  the  specific  volume,  P the  perssure,  S the  salinity  and 
V.  and  P0are  temperature  dependent  parameters  for  pure  water.  In  this 
per  the  thermal  expansion  data  of  Bradshaw  and  Schleicher  and  the  sound 
velocity  data  of  Wilson  have  been  incorporated  to  yield  a new  equation  of 
state  for  sea  water.  Values  of  the  adiabatic  gradient  calculated  with 
this  equation  for  S * 35*/*.  are  slightly  lower  at  2*C  than  those  reported 
by  Fofonoff,  for  example,  4%  at  400  bars  and  7%  at  1000  bars. 


UNIVERSITY  OF  CALIFORNIA,  SAN  DIEGO 
NARINE  PHYSICAL  LABORATORY  OF  THE 
SCRIPPS  INSTITUTION  OF  OCEANOGRAPHY 
SAN  DIEGO,  CALIFORNIA  921*2 


ALALYTIC  EQUATION  OF  STATE  FOR  SEA  HATER 


F.H.  Fisher,  Robert  Bruce  Williams  and  0*E.  Dial,  Jr. 


Sponsored  by 

Office  of  Naval  Research 
N000 14-63 *A-0200-6002 
NR  260-103 
and 

National  Science  Foundation 
DES- 70-00094 A04 


SIO  REFERENCE  75-30 


1 December  1975 


Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  U.S.  Government 


Document  cleared  for  public  release 
and  sale;  its  distribution  is  unlimited 


F.  N.  SPIESS,  DIRECTOR 
MARINE  PHYSICAL  LABORATORY 


MPL-U-99.1/67 


it 


Fisher,  Willi-t*?,  Dial 


PREFACE 

This  work  was  completed  in  19  T0  and  is  a sequel  to  the 
earlier  worst in  which  the  Tumlirz  . tuition  was  selected  to 
represent  the  properties  of  pure  water  as  well  as  sea  water. In 
this  paper  numerical  results  have  been  obtained  for  the  K|  and 
K2  paraweters.  Milli-ro's  group  has  developed  a more  accurate 
and  wore  complex  equation  of  state*  for  sra  wat*r.  This 
equation  has  been  used  by  soae  of  our  colleagues  at  Scripns  and 
has  :.een  found  useful  in  representing  properties  of  se»  water. 

F.  H.  Fisher,  December  1,  1975 

fF.  II.  Fisher  and  O.E.  Dial,  Jr.,  "Equation  of  State  of  Pure 
Mater  and  Sea  Mater,"  SIO  Reference  75-28,  l November  (1975). 

•F.  J.  Milltro,  R.  A.  Fine,  and  D.  P.  Mang,  "The  Equation  of 
State  of  Seawater,"  J.  Mar.  Res.,  32,  pp.  433-456,  (1974). 
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ABSTRACT 

Earlier  work  by  Fisher  and  Dial  on  pure  water  and  sea  water  showed  the 
effectiveness  of  the  Tumi irr  equation,  in  representing  the  precision 
pure  water  data  of  Kell  and  Mh.il ley  fro*  0 to  150*  up  to  1000  burs.  In 
addition  it  was  also  shown  that  a Modified  Tumlirz  equation  could  be 
used  to  represent  the  sea  water  data  of  Wilson  and  Bradley..  The 
equation  is 


V * Y'a,  -Kj  S ♦ 


P 


o 


V 

♦ K,  S ♦ 


I* 


where  V is  the  specific  volume,  I’  the  pressure,  S the  salinity  and 
V„,  X and  P0  arc  temperature  dependent  parameters  for  pure  water.  In 
this  paper  the  thermal  expansion  data  of  Bradshaw  and  Schleicher  and 
the  sound  velocity  data  of  Wilson  have  been  incorporated  to  yield  a new 
equation  of  state  for  sea  water.  Values  of  the  adiabatic  gradient 
calculated  with  this  equation  for  S = 35*/M  arc  slightly  lower  at  2®C 
than  those  reported  by  Fofonoff,  for  example,  4*.  at  400  bars  and 
7«  at  1000  bars. 


INTRODUCTION 

Fisher  and  Dial**0-  have  re-examined 
the  Tumlirz  equation  discussed  by  Eckart3  as 
an  equation  of  state  to  represent  pure  water 
and  sea  water.  Several  equations  were  examined 
and  the  analysis  included  data  from  several 
investigators  including  the  recent  precision  PV 
data  of  Kell  and  Nhallcy''  on  pure  water 
vrea  0 to  150*  and  up  to  1000  bars.  It  was 
found  that  the  Tumi  in  equation  for  specific 
volticr  V as  a function  of  pressure  P 


V » V,  ♦ p-Vl?  H) 

o 


provided  a fit  of  —8  ppm  from  0-100®  and 
about  IS  ppm  from  100-150®. 


In  attempting  to  apply  this  equation  to 
the  sea  water  data  of  Mi Ison  and  Bradley  5,6 
(they  also  used  the  Tumlirz  equation)  Fisher 
and  Dial1'2  found  an  erratic  dependence  of  P0, 
Va  and  X as  a function  of  salinity,  S as  shewn 
in  Figures  la,b  and  c.  This  was  also  true  of 
the  Newton  and  Kennedy  data7  . However,  when  X 
was  constrained  to  be  the  value  for  pure  water 
it  was  found  that  both  Pc  and  Vm  displayed  a 
linear  dependence  on  salinity  as  shown  in 
Figures  2a  and  b.  The  linear  dependence  of  the 
pressure  term  on  S had  been  suggested  by 
Tammann8  . 

Difficulties  were  encountered  in 
establishing  analytic  functions  for  the 
temperature  dependence  of  K j and  K2  solely  from 
the  Mi  Ison  and  Bradley  density  data,  partly 
because  the  data  for  various  salinities  were 
obtained  at  different  temperatures.  Also,  the 
results  for  thermal  expansion  coefficients  did 
not  agree  with  those  obtained  by  Bradshaw  and 
Schleicher  9 . 

In  this  paper  these  difficulties  have 
been  overcome  by  making  use  of  the  recent 
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theraal  expansion  data  of  Bradshaw  and 
Schleicher  * to  deteraine  K)  and  by  using 
Nil son's  sound  velocity  data10  to  deteraine 
Ka;  in  the  K 2 calculation  the  recent  specific 
heat  data  of  Bromley  et  al,11  are  also  ustu. 


Figure  la. 


O H>  to  M to 


Figure  lb. 


figure  lc. 

Figures  la,-  b,  and  c.  Plots  of  K , PQ  and 
for  best  (it  of  Tumlirz  equation  to  sea  water 
data  at  about  10*. 


The  equation  presented  here 
incorporates  the  aost  extensive  PV  data  on  sea 
water  froa  a single  laboratory  (see  Table 
11)5.1.7,12,13  as  wc)i  M tj,e  latest  theraal 
expansion  daca.  In  addition  the  equation 
provides  a high  precision  analytic 
representation  for  pure  water  when  the  salinity 
is  zero. 

Determination  of  Kj  and  K, 

The  full  equation  used  to  represent 
pure  water  froa  0 to  ISO*  is  shown  in  Table  I. 
Although  the  fora  of  equation  (1)  was 
determined  , froa  fits  to  the  Nilson  and  Bradley 
data,  difficulties  arose  in  trying  to  deteraine 
Kj  and  Kj  froa  the  volume  data  alone  since 
changes  in  Ki  could  be  compensated  for  by 
changes  in  K2.  However,  a determination  of  IC2 
can  be  made  at  atmospheric  pressure  using  sound 
velocity  data;  this  is  independent  of  the  PV 
data  as  w£ll  as  Kj*  The  isothermal 
compressibility  is 


3 ' - Htf),  ' ? 4 T (tt)p  S V 


where  V again  is  specific  volume,  c the  sound 
velocity,  T the  temperature  (K)  and  Cp  the 
specific  heat  at  constant  pressure.  Froa 
liquation  1 we  find  for 


, , 1 \ 

® V 2 * 

(P0  ♦ K2  S ♦ P)Z 
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Figur e 2a. 


Figures  2a  and  b.  Plots  of  P„  and  V for  best 
fit  of  Tumlirz  equation  to  sea  water  data  at 
about  10 * when  is  constrained  to  be  value 
for  pure  water. 


and,  therefore,  at  1 bar 


*2  S 1 -<Po  * *> 


1 X 

|v2/c2*T 

cV' 

l3T< 

1 /cr 
v,s 

(4) 


The  temperature  dependence  of  K was  determined 
by  using  data  for  Cp  from  Bromley 9 , sound 
velocity  from  Wilson10,  specific  volume  from 
Wilson  and  Bradley  and  (SV/  3T)  from  Bradshaw 
and  Schleicher  9 . 


Next,  to  determine  Kj  the  thermal 
expansion  data  of  Bradshaw  and  Schleicher  were 
used.  The  general  form  of  K|  and  Koas 
functions  of  temperature  were  determined  from 
the  Wilson  and  Bradley  1’V  data  to  be 


K ,*  B ♦ B,  T ♦ B,  T2 
I O 1 2 


and 


K2  = no  * h T 


(5) 


(6) 


By  differentiating  equation  (1)  with 
respect  to  temperature  the  only  unknown  term  as 
a function  of  temperature  is 


= B. 


2B,  T 


(7) 


Tli  is  was  evaluated  by  using  the 
Bradshaw  and  Schleicher  data  fori  3V/3T)pa<;. 
After  this  the  Wilson  and  Bradley  data  was  used 
to  find  the  best  values  for  Bu. 

RIISUI.TS 

The  temperature  dependence  of  Va, 
P0  and  X determined  from  the  pure  water  data 
are  shown  in  Table  I for  both  the  entire  0-150° 
range  of  the  Kell  and  hhallcy  data  and  in  a 
truncated  form  that  is  limited  to  the  same 
temperature  range  of  sea  water  data.  The 
temperature  dependence  of  both  Kj  and  K2  arc 
also  shown  in  Table  I.  The  coaplctc  equation 
for  pure  water  exhibits  a density  maximum  at 
4.00®. 

Table  III  shows  the  fit  of  our  equation 
in  Table  I to  the  Wilson  and  Bradley  data. 
Although  the  overall  fit  is  within  the 
precision  claimed  by  Wilson  and  Bradley  it  is 
•ecn  that  the  low  temperature  data  shews  poorer 
fits  than  the  room  temperature  data. 

Table  IV  summarizes  results  of  various 
investigators  for  the  thermal  expansion 
coefficient  of  sea  water  shown  in  the  Wilson 
and  Bradley  technical  report  along  with  our 
results. 

Table  V in  an  analogous  manner 
summarizes  the  specific  volume  results. 

Similarly,  Table  VI  summarizes  the 
compressibility  results. 

In  Table  VII  we  compare  our  results  for 
thermal  expansion  ( 3 V/  3 T)j>  5 against  those  of 
B5S  as  a function  of  temperature,  salinity  and 
pressure.  At  0"  we  find  the  greatest 
disagreement  with  their  results, 
~10'.  at  1 bar  and  ~3'#  at  1000  bars. 

In  Table  VIII  we  do  see.  however,  that 
changes  in  specific  voIimcs  from  -2°to*2° 
calculated  by  our  equation  compare  very 
favorably  with  the  results  of  BtiS.  In  this 
case  we  arc  comparing  values  from  our  equation 
against  their  raw  data. 
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Table  11 

Comparison  of  Density  Data  at  elevated  Pressures  for 
Pure  Mater  and  Sea  Mater  in  Teaperaturc  Range  O'  - 40* , 


Same 

Tenp/Pressure 
kance/ Isotherm  Precision 

Number  of  Data  Points  by 

Total 

Data 

Points 

Tenp 

Pressure 

Salinity 

Pure 

Mater 

EKaan 

0 - 20*/  600  Bar 

5 

3 

2 

1 

45 

Newton- Kennedy 

0 - 2571300  Bar/70  ppm 

6 

13 

3 

1 

312 

Milson-BraJlcv 

0 - 40*/  965  Bar/ 20  ppm 

9 

15 

s 

l 

795 

Kell-Mhalley 

0 - 40°/ 1026  Bar/5  ppm 

8 

27 

0 

1 

216 

Amagat 

0 - 4071000  Bar 

7 

10 

0 

1 

70 

3 

Z 


Table  111 


Standard  deviations  in  parts  per  million  of  the  fit  of  the  FMD 
equation  to  raw  specific  voltMc  data  of  Milson  and  Bradley  vs 
tenperature  and  salinity.  (Temperature  cited  is  nominal  since 
data  were  taken  at  varying  temperatures  near  cited  one.)  The 
average  of  the  standard  deviations  is  10"  ppm  for  all  the  data. 


•MM 


^ C/..) 

10.221 

20.720 

30.881 

35*5t>S 

40.3*0 

- 0° 

208 

175 

250 

264 

156 

- 5° 

1*3 

129 

212 

1S7 

236 

-10° 

128 

5! 

16" 

!"2 

45 

-15r 

102 

58 

94 

134 

55 

—20° 

42 

76 

44 

140 

37 

~2S3 

17 

60 

85 

95 

62 

-30* 

17 

106 

95 

47 

104 

-35° 

22 

11" 

55 

45 

98 

~40* 

19 

154 

54 

100 

161 

t ’) 

‘ ? 
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Tabic  IV 

Thermal  Expansion  of  Sea  Kater  (S  * 35*/.. ) 
Comparison  Tables 


u 

P 

Bars 

Hydrographic 

Tables 

x IO3 

Eckart 
x 105 

Crease 
x 105 

SOI 
x 105 

FKD 
x 103 

0* 

1 

5 2 

8.0 

5.6 

7.8 

5.76 

200 

10.5 

13.4 

10.6 

15.2 

10.9 

400 

15.4 

18,2 

15.0 

18.0 

15.4 

600 

19.8 

22.5 

17.1 

22.2 

19.5 

. 

BOO 

25.2 

26.6 

22.6 

26.1 

22.6 

1000 

26.5 

29.5 

25.6 

29.5 

25.5 

io® 

! 

16#  7 

16.2 

16.7 

16.5 

16.6 

200 

20.2 

19.9 

20.1 

20.0 

20.1 

: 

400 

23.3 

23  1 

23.2 

23.5 

25.2 

600 

26.1 

26.0 

26.1 

26.5 

25.9 

800 

28.6 

28  8 

28.6 

29.0 

28.3 

1000 

30.9 

50.“ 

30.  S 

51.5 

30.5 

20® 

I 

25.7 

23.8 

25.7 

24.1 

25.6 

200 

27.8 

26.2 

27.7 

26.5 

27.8 

400 

29.8 

28.2 

29.7 

28.6 

29-7 

600 

51.7 

30.0 

50.6 

31-5 

800 

55.5 

51.7 

32.3 

55. 1 

1000 

35.5 

32.9 

55. 8 

34-6 

30* 

I 

35  5 

51.5 

55  • 4 

31-7 

55.1 

200 

54-  7 

52. 5 

54.6 

55  0 

54.5 

400 

36. 1 

55.6 

55  7 

34  2 

35-6 

600 

37.6 

•4,4 

55  2 

56.5 

800 

59.5 

55.2 

56.  I 

57  5 

1000 

41.4 

35.8 

36-9 

38.5 

40* 

I 

41.0 

39.4 

59.6 

40.5 

200 

41.8 

59.5 

40.9 

40.7 

400 

43.0 

59.5 

40.2 

41.1 

600 

44.6 

3°  5 

40-5 

41.4 

800 

46.9 

39.  .> 

40  6 

41.7 

1000 

49.6 

39.5 

40.7 

42.1 
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Tabic  VI 


Compressibility  of  Sea  Mater,  S * 55°/*.  (Bars**) 
Co^urison  Tables 


n 

P 

8ars 

Hydrographic 

Tables 

* 10b 

fc'ckart 
x in6 

Crease 
x 10b_ 

N0L 
x 106 

FND 
x 106 

0° 

i 

46.2 

46.5 

46.4 

46-7 

46.3 

200 

44.0 

44  « 

44.0 

44.1 

45.9 

400 

41  7 

41.7 

41.7 

41.7 

41.7 

600 

39.7 

59-5 

39.6 

39.5 

39.6 

800 

37. S 

57.  S 

37.5 

37.7 

1000 

36.2 

55*7 

55.9 

35.7 

36;0 

10 1 

1 

44.1 

44.1 

•54 . 1 

44.3 

44.0 

200 

42.0 

41.8 

41.9 

42.0 

41.9 

400 

39.9 

39.7 

39.9 

59.8 

39.9 

600 

38.0 

37.  8 

5S.1 

37.8 

38.1 

500 

36-4 

56.0 

36.5 

36.0 

56.4 

1000 

34.9 

51.5 

54.6 

54,5 

34.8 

:o’ 

1 

42.6 

42  5 

42.7 

42.7 

42.6 

200 

40.6 

40.5 

40.6 

40.5 

40.6 

400 

58  7 

58*5 

3S.7 

38.5 

38.8 

600 

56  9 

36.7 

56*7 

37.0 

•00 

5575 

35.0 

35.0 

35.4 

1000 

55.9 

35.4 

55.4 

33.9 

30’ 

1 

41  S 

41.6 

41.9 

41.7 

41.7 

200 

39.  S 

59.6 

59.8 

39.7 

39.8 

400 

3*  .9 

57,  S 

38.0 

37.8 

38.1 

609 

56. 1 

M*  I 

36.0 

36.4 

MO 

54.4 

54.4 

34.4 

34.8 

1000 

32.9 

52.9 

52.8 

33.4 

40 1 

1 

41  4 

41.5 

41.2 

41.4 

200 

59.5 

39.5 

39.3 

39.5 

400 

37.2 

57  5 

37,4 

37.7 

600 

55*5 

35.8 

35.7 

36.1 

- 

800 

55.4 

34.5 

34.1 

34.5 

1000 

31.6 

52  8 

32.6 

33.1 
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Table  VII 


Coeparison  of  Values  of  in  parentheses  calculated  fro*  FWD  Equation  with  Those  of 

' c T/P,S 

Bradshaw  and  Schleicher  Co^Kited  free  Their  Comprehensive  Foraula. 


S « 30.50*/* 


(1S4)  154 
(230)  229 
(284)  284 


S * 35.00*/, 


P.bar?^ 

0 

10 

20 

30 

1 

(56) 

52 

(162) 

162 

(250) 

251 

(32^) 

327 

500 

(166) 

166 

(234) 

234 

(293) 

295 

(346) 

347 

1000 

(238) 

244 

(285) 

286 

(325) 

325 

(361)  _ 

363 

S * 39.50*/, 


XQC 

P.bar^N. 

0 

10 

20 

30 

1 

(68) 

65 

(169) 

170 

(2S4) 

256 

(329) 

329 

500 

(172) 

174 

(238) 

239 

(295) 

296 

(347) 

348 

1000 

1241) 

248 

(287) 

289 

(325) 

326 

(361) 

363 

Fisher,  Williams,  Dial 


In  Table  IX  we  show  a comparison  of 
sound  speeds  calculated  from  our  equation  with 
those  reported  by  Del  Grosso14  along  with  Del 
Grosso's  comparison  of  his  results  with  others. 
Figure  3 shows  a plot  of  these  differences 
which  indicate  that  sound  speed  calculated  from 
our  equation  is  in  general  agreement  with  all 
data  except  the  latest  from  Del  Grosso. 

DISCUSSION 

Of  special  interest  to  physical 
oceanographers  for  stability  calculations  is 


the  adiabatic  gradient.  Values  of  the 
adiabatic  gradient  calculated  equation 
for  S = 35  were  slightly  lower  than  those 
reported  by  Fofonoff  ls,  for  example,  4Vat  2* 
and  £00  bars  and  7*  at  2s  and  1,000  bars.  In 
addi.  ion  we  also  find  that  the  density  maximum 
crosses  the  freezing  point  at  S*  22*/.» 
inst  >ad  of  the  usually  cited  value  of 
S 5 25V... 

Since  Li's  equation  did  not  include 
ah  expression  for  the  density  of  water  at 
atmospheric  pressure  we  have  not  calculated 
thermal  expansion  coefficients  and  other 


Table  VI 11 


Calculated  values  from  Fh'D  equation  vs  BfiS  observed  values 
for  the  change  in  the  specific  volume  of  sea  water  from  -2* 
to +2*.  Fh'D  values  are  listed  in  parentheses:  Units  arc  in 


cc/gm. 


P bars 

Salinity  •/•• 

201.3 

. 25.0041.  _ . 

401.2 

35.005 

601.0 

35.004 

800.0 
. „ 35.002 

mMi 

- - - 

TJC 

(-260) 

(-352) 

(-415) 

(-467) 

-2 

(-101) 

-277.1 

*o5(>.9 

-424.3 

-480.5 

-1 

-9T.5 

(0) 

(0) 

(0) 

(0) 

(0) 

0 

0 

0 

0 

0 

0 

(230) 

(310) 

(370) 

(4371 

(486) 

-2 

225 

310 

383 

445 

498 

Table  IX 


Comparison  of  Speed  of  Sound  Values  calculated  from  the  FWD 
liquation  of  State  with  those  reported  by  Del  Grosso  for  NRI. 
and  .SOI.  data. 


T 

Nil  son  (1960) 
N0I.-2 

Del  Grosso  (1970) 

FWD 

.141 

1449.05 

1449.61 

1449.86 

10.000 

1490.40 

1489.75 

1489.98 

19.036 

1521.96 

1521.13 

1521.90 

29.970 

1545. "7 

1545.42 

1546.39 

40.027 

1564.77 

1563.11 

1564.37 

properties  from  his  equation.  We  did  check  7. 
som  specific  volumes  calculated  using  his 
-.equation  and  they  were  within  166  ppm  of  ours  8. 
at  30\ 
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